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In this note we discuss emergence of geometrical scaling (firstly proposed for deep 
inelastic collisions) in pp scattering at the LHC and in heavy ion collisions at RHIC. 
After discussing general properties of geometrical scaling (GS) we focus on simple 
signatures of GS, namely on ratios of the pt values at which multiplicity spectra 
for different energies are equal. 

1 Geometrical Scaling in Deep Inelastic Scattering 

One of the outstanding problems of the evolution of parton densities is the rapid increase 
of the number of gluons at low Bjorken x lj. Such growth has to be tamed at some 
point. The scale at which this happens is called saturation scale Q s and it depends on 
the Bjorken x. The explicit form of the saturation scale follows from the fact that Q 2 (x) 
is related to the gluon distribution in the proton at low x 



where Qo ~ 1 GeV and xo ~ 10 -3 are free parameters whose precise values can be 
extracted from fits to the HERA data. Power A is known to be of the order A ~ 0.2-^0.3. 

Based on the saturation idea Golec-Biernat and Wiisthoff proposed a model [2] where 
the reduced cross-section cr 7 * p ~ F2/Q 2 is written as: 



where tp is an amplitude for photon dissociation into a qq pair of size r, which is explicitly 
known, and a^p is a model dipole-proton cross-section. 



Q 2 s {x) = Ql (x/x Q y x 
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There are two important features of Gdp that are incorporated in the GBW model. 
Firstly, instead of taking a^p to depend independently on r and W (7*p CMS energy) 
it is assumed that <JdP depends on Bjorken x 

x = Q 2 /(Q 2 + W 2 - M 2 p ) (3) 

where M p stands for the proton mass. Furthermore, in the GBW model o^p depends on 
x through the saturation scale Q s ([TJ, namely: 

<J d p{r 2 ,x) = a dP {r 2 Q 2 (x)). (4) 

Given the fact that for massless quarks distribution amplitude tp is a function of a product 
of r and Q: 

mr,Q 2 )\ 2 = Q 2 mrQ)\ 2 (5) 

it follows that 

cr rp = function(Q 2 / ( 9s( a; )) (6) 

a phenomenon known as geometrical scaling [3 . 

Let us finish this section by recalling that the dipole-proton cross-section a d p is in 
the GBW model related to the unintegrated gluon distribution ip(x.k^) [4j[2]: 

4tt 2 fdkl 



Vd P {x, r) = — J -p- (1 - J (rfc T )) a s ip(x, fc T ) (7) 

where 



Q 

xG(x,Q 2 ) = j ' dk\ip{x,k^). (8) 
The fact that a dp cross-section depends only on the combination r 2 Q 2 (x) implies that 

l p(x,k 2 ) = v (k 2 /Q 2 (x)) (9) 
and strong coupling constant a s should be frozen. 

2 Geometrical Scaling in pp Collisions 

For pp collisons for low and moderate pt one uses Gribov-Levin-Ryskin formula [5] 

j ^ = 2 f d 2 k T a s (p 1 (xi,k 2 : )(p2(x2, (fc-p)|). (10) 
drjd 2 p T p^J 
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Here Xi t 2 are gluon momenta fractions 

Xl2 = e ±v p T /W with W = -/a (11) 

needed to produce a gluon of transverse momentum and rapidity y (or pseudorapidity 
rj) . Unfortunately formula ( 10 ) has been proven jS] only for the scattering of a dilute 
system on a dense one (i.e. x\ << x 2 or 12 « x\). Despite that, one is often forced to 
use ( |10| ) in a region of X\ ~ where numerical studies show that it still works reasonably 
well (see e.g. Ref.[7]). 

If unintegrated gluon densities scale according to Eq. Q than also dN/ dr]d 2 pT should 
scale (provided we freeze a s and neglect possible energy dependence of constant C). 
Therefore geometrical scaling for the multiplicity distribution in pp collisions [H [10] 
states that particle spectra depend on one scaling variable 

t = P 2 t /QHpt,W) (12) 



where Q 2 (pt, W) is the saturation scale at x\ ~ x 2 ( 11 h 

Q 2 s (PT,W) = Q 2 (p T /(WxlO- 3 )y X (13) 

where we have neglected rapidity dependence of 2:1,2 • Factor 10 -3 corresponds to the 
choice of the energy scale (arbitrary at this moment xq in Eq.([T])). Hence 

dof dNch 



N ch (W, PT ) 



di]d 2 pT 



IsHr) (14) 



with Qq ~ 1 GeV. Here F(t) is a universal function of r. 

In order to examine the quality of geometrical scaling in pp collisions in Ref.|ll| we 
have considered ratios R\Vi/w 2 

p , v dcf N ch (W uPT ) 
N ch (W 2 ,PT) 

Here we shall discuss another way of establishing geometrical scaling, at least quali- 
tatively. Note that if at two different energies W\ and W2 multiplicity distributions are 
equal 

N&iWupQ) = N ch (W 2 , P (2) ) (16) 



then this means that they correspond to the same value of variable t (12). As a conse- 
quence 

^ )2 (^V^) A =P? )2 (^ ) /W 2 ) A (17) 
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for constant A. Equation (171 implies 



OPT def (1), (2) 

lD w 1 /w 2 ~ IPt 



(Wi/W 2 ) 



(18) 



Ratios &w i , Ws for pp non-single diffractive spectra measured by the CMS [3] col- 
laboration at the LHC are plotted in the left panel of Fig. [I] together with the straight 



horizontal lines corresponding to the r.h.s. of Eq.(18) for A = 0.27. We see approximate 
constancy of &w 1 /w 2 over the wide range of iV c h 



A small rise of S^^ w ^ with decreasing 



iV c h corresponds to the residual px-dependence of the exponent A |10j 



S ratios 7/0.9 and 7/2.36 
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Figure 1: Plots of S^/wv Left panel: CMS pp spectra for W x = 7 TeV and W 2 = 
0.9 TeV (blue full circles) and 2.36 TeV (red triangles) in function of log(iV c h). Right 
panel: PHOBOS Au-Au spectra for Wi = 200 GeV and W 2 = 62.4 GeV in for two 
centrality classes: — 5% (red triangles) and 45 — 50% (blue full circles) in function of 



\og(N cil ). Straight lines correspond to the r.h.s. of Eq. (18) for A = 0.27 



We would like to advertise this method as a simple way of looking for GS in the 
Pt spectra in various reactions. An obvious advantage is that it is very easy to do: it 
requires only to interpolate one spectrum (here we have done this for W\ = 7 TeV). An 
obvious disadvantage consists in the fact that it is very difficult to attribute sensible error 
to the ratios S^, W2 , so for quantitative purposes it is better to consider ratios R\y 1 /w 2 - 
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3 Onset of Geometrical Scaling in Heavy Ion Collisions 

We have already briefly discussed the onset of GS in heavy ion collisions in Ref. [TT] . 
This work is still in progress. If GS in HI collisions holds then, analogously to equation 



(18), we can form ratios of pr's corresponding to the same multiplicity (in the same 
centrality class c) at two different energies W\p,. In what follows we will neglect possible 
dependence of the saturation scale upon participant density [12] (for review see Ref. [13]). 

Let us concentrate on the PHOBOS data for Au-Au [HI [15]. PHOBOS data cover 
asymmetric rapidity region 0.2 < rj < 1.4 at two RHIC energies 64.2 and 200 GeV. The 
highest centrality bin for PHPBOS is only 45-50 %. In the right panel of Fig. [I] wc 
plot ratio S'200/62 4 ^ or two centrality classes: 0-6 % and 45-50 %, together with a line 
corresponding to A = 0.27. From Fig.[l]we might conclude that geometrical scaling works 
well also in heavy ion collisions. More detailed studies show, however, that dependence 
of the saturation scale on the participant densities may be crucial for proper description 
of the data. 



4 Conclusions 

In this note we have proposed a simple procedure to look for geometrical scaling in the 
Pt spectra, namely to construct ratios of transverse momenta corresponding to the same 
multiplicity. Based on this method one can qualitatively see GS in the CMS spectra [9] 
at the LHC energies and in HI collisions at RHIC fTU [T5] . 

Many aspects of GS require further studies. Firstly, new data at higher energies 
(to come) have to be examined. Secondly, more detailed analysis including identified 
particles and rapidity dependence has to be performed. On theoretical side the universal 
shape F(t) has to be found and its connection to the unintegrated gluon distribution 
has to be studied. That will finally lead to perhaps the most difficult part, namely to 
the breaking of GS in pp. 

Even more work is needed to understand GS quantitatively for HI collisions. Here 
the main question arises why hydrodynamic evolution preserves GS over a time of a few 
fermi until the final state particles are produced. Moreover, the behavior with energy of 
transverse participant densities for different centralities and for different rapidity regions 
require further studies. And finally A dependence of the saturation scale has to be 
confirmed. Here the new LHC data for different energies and different nuclei will open a 
new chapter, since the LHC energies are order of magnitude higher than those of RHIC. 
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